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Abstract 

The geometric entropy arising from partitioning space in a fluid 
"field theory" is shown to be linearly proportional to the area of an 
excluded region. The coefficient of proportionality is related to sur- 
face tension by a thermodynamic argument. Good agreement with 
experimental data is obtained for a number of fluids. The calculation 
employs a density matrix formalism developed previously for studying 
the origin of black hole entropy. This approach may lead to a prac- 
tical new technique for the evaluation of thermodynamic quantities 
with important entropic components. 
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1. Prolegomena 

The ground state of a quantum field theory can be described by means of a 
density matrix. If the field degrees of freedom inside a specified volume are traced 
over, the result is a reduced density matrix, p out , which depends only on the degrees 
of freedom external to the excluded volume. Associated with this reduced system is 
a "geometric entropy", S out /k = —Tr(p out \np out ), which quantifies the information 
lost by the partitioning of space. The geometric entropy is an expression of the fact 
that p out lacks the information contained in correlations between the excluded interior 
volume and the rest of the system. 

This density matrix formalism was recently utilized in an attempt to provide a 
simple explanation of the classic result M that the entropy of a black hole is linearly 
proportional to its surface area. A spherical volume of a space containing a free scalar 
field was excluded, and the resultant geometric entropy was then determined. It was 
found' 2 ' 3 ' that the geometric entropy of this "black hole" is, in fact, proportional to 
the area of the excluded region, rather than to its volume. This result seems initially 
to be rather mysterious, since entropy (like free energy) is generally an extensive 
quantity. However, in several important physical situations, major contributions to 
the free energy of a system are, in fact, proportional to its area. These "area law" con- 
tributions are primarily responsible for liquid surface tension and for "hydrophobic" 
effects central to protein folding. 

Although the density matrix formalism was developed to further our understand- 
ing of black holes, it does not directly utilize either classical general relativity or 
quantum gravity. It will be shown here that geometric entropy can be used to exhibit 
a point of commonality between liquid surface tension and the aforementioned results, 
which are germane to the study of black holes. The geometric entropy arising from 
partitioning space in an empirical fluid "field theory" is calculated. Remarkably, the 
geometric entropy is found to be linearly proportional to the area of the excluded 
volume, as in the above case' 2 ' 3 ' of a free scalar field. A thermodynamic argument re- 
lates the coefficient of proportionality to liquid surface tension. Good agreement with 
experiment is obtained for a number of liquids, suggesting that the density-matrix 
formalism may lead to a simple and accurate new way to evaluate thermodynamic 
quantities with important entropic components. Further anticipated developments 
are mentioned in the closing remarks. 
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2. A field theory of fluids 



The applicability of the density-matrix formalism to fluids is founded upon the 
observation that a liquid at finite temperature can be considered, for some purposes, 
as the vacuum state of a field theory This is a natural step to take, since the 
structure and properties of a fluid at equilibrium can be described by its molecular 
distributions' 4 ' 5 ^ . These distributions can, in turn, be used to define the Green's 
functions of a field theory, e.g., 

<P(r')> = Po 

(p(r')p(r' + r)> = pg[l + h(r)\ + p 5(r) (1) 

where the right-hand sides of Eqs. (1) are independent of r' for a homogeneous fluid. 
These Green's functions (• • •) are vacuum (i.e., ground state) expectation values of 
fluid density fields (p(r)}, calculated by means of a functional integral measure defined 
below. The complete set of Green's functions, which can be taken as the definition of a 
field theory, can likewise be expressed in terms of higher-order molecular distributions. 
The average fluid number density po an d pair correlation function h(r) are physically 
measurable, and can be extracted from experimental data. The first two Green's 
functions are thus known quantities. 

The idea of representing a fluid at finite temperature by the vacuum state of a 
field theory is quite novel, and therefore deserves careful explanation. Relations such 
as Eqs. (1) are frequently found in the classical treatments of fluids' 4,5 ]. However, in 
the "classical" treatment, the averages (• • •) are meant to be taken over canonical or 
grand canonical ensembles. Thus, although the definitions Eqs. (1) are identical to 
equations commonly found in textbooks, their interpretation is quite different. 

In the classical formalism, one usually begins by taking averages over a canonical 
ensemble of N particles which interact via a potential function V(ri, r 2 , • • • , tn). Then 
the canonical ensemble averages (• • •) canonical of densities and their correlations are 
defined by relations like 

(p( r ))canomcal = Z ^ X J *1 J dr 2 ■ ■ ■ J dr N exp(-/3F) X p dassica l (r) 

where 
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N 

Pclassical(r) = ^ 5(r - n) 
i=l 

The integration is over the coordinates of all particles 1 through N, (3 = l/fcT is the 
inverse temperature, and Z^v is the partition function. Grand canonical averages can 
then be calculated in the usual fashion to yield the familiar molecular distributions. 

One can thus generate an infinite set of grand canonical expectation values 

(p(ri)p(r 2 )---p(r M )) gc 

for all products of the quantity p c iassicai( r ) ■ It is then that a novel step is taken. Each 
grand canonical average (• • •) is identified with the corresponding Green's function of 
a certain field theory. Because the complete set of grand canonical expectation values 
is known in principle, all of the Green's functions of this field theory are also known. 
The full set of Green's functions for a field theory can then be taken as a definition of 
that field theory. One important point must be emphasized, however: field theoretic 
Green's functions are, by definition, expectation values of certain operators taken over 
the vacuum state. Thus, if the operators and field-theoretic ground state are chosen so 
as to reproduce the hierarchy of grand-canonical expectation values discussed above, 
the ground state of a field theory provides a description of a fluid that is precisely 
equivalent to knowledge of the complete set of molecular distributions. 

In practice, of course, neither field theories nor grand-canonical fluid models are 
generally solvable. Thus, one must resort to judicious approximation. As the present 
study is primarily intended to illustrate the potential benefits of a very new approach, 
the field theory used is a simple one, a generalized free field theory chosen to reproduce 
only the first two Green's functions Eqs. (1). 

The simplest field theory whose Green's functions reproduce Eqs. (1) has a 
ground-state density matrix proportional to exp[— ^(S{p}+S{p'})], where S is quadratic 
in the field variables {p(r)}: 

PoS{p} = l -JdvJdv'[p{v)-p Q }W{v-v')W)-p,} 
W{y-t') = 5(r - r') - p c(r - r') (2) 
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and the function c(r — r') remains to be determined. As an additional convenience, 
the field variables {p(r)} are allowed to range from negative to positive infinity. The 
expectation values in Eqs. (1) are calculated via the functional integral 



(O) = Z- 1 Jo{p}eMS{p}]Vp 
S = jeMS{p}]Vp 



(3) 



The function c(r) is then required via Eq. (1) to satisfy 




(4) 



which is the well-known Ornstein-Zernike relation' 6 ' linking the direct correlation 
function c(r) with the pair correlation h(r). The physical meaning of the function 
c(r) is thus manifest. The success of the Ornstein-Zernike approach derives in part 
from the short-ranged character of c(r), suggesting that Eq. (2) is, indeed, a sensible 
starting point. 

For this choice of c(r), Eqs. (2) and (3) reproduce Eqs. (1). The constraints 
posed by Eqs. (1) alone are, however, insufficient to define a field theory uniquely. 
The generalized free field theory implicit in Eqs. (2) and (3) is only the simplest 
solution of Eqs. (1), and therefore provides an approximate description of a real 
liquid. In the theory defined here, all higher-order correlations are calculable in 
terms of products of those in Eqs. (1). The fluid theory defined by Eqs. (2) possesses 
the property that higher-order molecular distributions factorize, which is known' 5 ! 
to be correct for large separations. There is no barrier of principle to performing 
a more elaborate calculation with a better model of a fluid; however, the labor of 
calculation would be larger. (This point is discussed in further detail below). The 
integration over unphysical negative values of the fluid density fields {p(r)} is also 
permissible, for S{p} is positive definite and is sharply peaked about p(r) = p > 0, 
provided that the Fourier transform W(k) is positive (as it is for a fluid with finite 
compressibility). Although negative, "virtual", values of the density appear formally 
in the functional measure Eq. (3), negative densities do not explicitly enter into the 
physically observable molecular distributions, which are the Green's functions of the 
theory. The standard' 5 ' density functional approach can be recovered by adding the 
usual potential term J v(r)p(r)dr to S{p}, and then evaluating S as a functional of 
the average densities p(r) = (p(r)) v . The density functional fl{p} of the bulk fluid is 
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then simply f2{po} — kT x ln[S{p}], which is the "effective action" of the field theory 
Eq. (3). For any S{p}, fi{p} calculated in this fashion is convex^ for a uniform 
system and can have no more than one minimum, so no "Maxwell construction" is 
needed. 

In what follows, the fluid density po and pair correlation function h(r) are taken 
from experimental results. No attempt is made to derive either of these quantities 
from first principles. One final simplification is therefore employed in the sequel. 
Experimental techniques (such as x-ray scattering) typically measure only a spher- 
ical average h(\r\) of the pair correlation. Although the above formalism can be 
easily applied to test a nonspherically symmetric model, using real data requires 
the assumption that h(r) = h(\r\) and c(r) = c(|r|). This assumption is, however, 
commonplace' 12 ' . 

3. Relating geometric entropy to surface tension 

The geometric entropy can be related to fluid surface tension by a simple argu- 
ment. A well-known theorem of statistical mechanics' 8 ^ states that the probability of 
observing a fluctuation in a system at temperature T is given by exp(— W/kT), where 
W is the reversible work required to produce the fluctuated configuration through the 
application of a constraint, and k is Boltzmann's constant. Thus the reversible work 
required to produce a cavity of low density in a fluid can be found from the proba- 
bility that such a cavity occurs via a fluctuation. Consider a subsystem of the fluid 
contained within an imaginary sphere of radius R. The density matrix p in for this 
subsystem can be constructed from that given above for the full volume by integrating 
out the field variables (p(r)} outside the sphere. This density matrix is a weighted 
sum of the states internal to the sphere, given that the exterior region is unobserved. 
It is shown below that, ignoring volume work, the relative probability of occurrence of 
any state which approximates a low-density cavity is essentially exp(—S in /k), where 
Sin/k = — Tr(pi n In pin) is equal to the geometric entropy S out /k of the subsystem. 
The macroscopic surface tension 7 for a fluid at a fixed temperature T is then the large 
R limit of TS in (AiiR 2 )~ l . Since the Green's functions used to specify the fluid field 
theory depend implicitly upon the liquid temperature, geometric entropy surface ten- 
sion varies nonlinearly with temperature, and thus contains contributions from both 
surface excess entropy and enthalpy. 

The useful correspondence between this surface tension calculation and the afore- 
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mentioned "black hole" entropy result arises upon construction of the complemen- 
tary density matrix p out . In this complementary case, only the fields (p(r)} inside 
the sphere are integrated out to yield p out . Then it can be shown' 3 ] that S in /k = 
Sout/k = —Tr(p out In p out ), where S out /k is the geometric entropy of the subsystem. 
The density-matrix calculations' 2 ' 3 ^ relevant to black hole entropy involve the devel- 
opment of techniques for calculating S out ; by the reasoning given here, they are also 
useful for extracting liquid surface tension. 

The density matrix calculation proceeds as follows. When the field variables {p(r)} 
outside the sphere are integrated out, the result ' 2 1 is a reduced density matrix 

p in {q,q'} = [det(M/n)^eM-l(QMq + q'Mq') - ±(g - q')N(q - q')} (5) 

where {q(r) = [p(r) — po]} inside the sphere. The matrix M is the inverse of W/ (2p ), 
taken over this interior domain, while N = W/ (2p ) — M. Since this domain is finite, 
M and N are to be considered sums over discrete values of wavenumber. 

The reduced density matrix p in is next expressed in a new basis {x} in which the 
matrix A5 = M~?NM~2 is diagonal, with eigenvalues A. The coordinates {x} are 
normalized by x 2 = qMq, so that the determinant prefactor in Eq. (5) is eliminated. 
Then' 2 ' 3 ! for each mode [with eigenvalue A = 4£(1 — £) -2 , < £ < 1]: 

00 

Pin,\{x, X'} = Pn^n(x)^ n (x') (6) 

?1=0 

* n (x) = (2 n n! x^J^c)-^H n (^x) exp(-ax 2 /2) (7) 

where H n (z) is an Hermite polynomial, a = (1 + £)/(l — £,), and p n — (1 — ^)^ ra . The 
reduced density matrix is a product of those for each mode A: 

p in {x, X'} = Y[ Pin,\{x, X'} (8) 
A 

Each mode Pi n ,\ behaves like a thermal density matrix for an harmonic oscillator, 
specified by a frequency a and effective temperature T e ff = a/ln(l/£). 

The probability that a mode A of the spherical subsystem is found in state n is 
thus p n . Associated with this probability distribution are two useful quantities. These 
are n, which is the average state number: 
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n=Y,nxp n = a(l-0 (9) 

?1=0 

and S in , the "geometric entropy" , in terms of which the density of states is exp(—S in / k): 

oo 

Sin/k = -Y / Pnlnp n = - ln(l - - nine (10) 

n=0 

The probability of finding the mode in any state n > N is p any {N) = £ N . Thus, 
for n > n, 

Pn<Pn = exp(-S in /k) < p any {n) (11) 

For the systems of interest here, £ is generally quite small, so the second inequality 
is nearly an equality: 

exp(-S in /k) = Pany(n) (12) 

In the classical methods of calculation, the probability that the system attains 
a given density is calculated. The course taken here is, however, conceptually quite 
different. The quantity of interest here is the relative probability to find the subsystem 
in certain states, those whose characteristic density is relatively low. For a given 
state, there is a finite probability that any value of the density will be attained. 
The expectation value of an operator can, however, be determined for such a state. 
Thus, the expectation value of x 2 ~ (^{p} — £l{po})/kT in state n is given by 
(n + \)/a. The probability distribution [\I/ n (a;)] 2 for n > typically is maximized 
when x = x pea k = ±[(2n + l)/a]^, corresponding to values of the density that are 
above and below the bulk value. A state whose density is expected to differ from 
the bulk value by at least the inverse volume of the subsystem (so that at least one 
molecule is absent) must therefore have n> n. By Eq. (12), the probability of finding 
the system in such a state is well approximated by exp(— Si n /k). 

Where, then, are the approximations made, and how does one deal with more 
complex situations? Let us review the logic of the calculation. One begins with a 
field theory in its ground state. The question of interest is the probability that a given 
fluctuation occurs inside an imaginary sphere of radius R, regardless of what occurs 
outside that imaginary sphere. Thus, the degrees of freedom outside the imaginary 
sphere are integrated out. What remains is a field theory that is defined in terms 
of variables internal to the sphere. This field theory gives the relative probability 
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of a given internal configuration, assuming that the region external to the sphere is 
unobserved. In this regard, there is a point of contrast between the fluid calculation 
and the black hole entropy result. In the latter case, degrees of freedom inside the 
black hole are considered unobservable on fundamental grounds, rather than simply 
being unobserved. Nevertheless, the correct formal procedure is the same in the two 
cases-the unobserved variables are integrated out. 

It is important to realize that, at least in principle, no approximations need to be 
made to reach this stage of the calculation. Instead of the quadratic form given in 
Eq. (2), one could just as well begin with an arbitrarily complex density functional 
S{p}, corresponding to an arbitrary set of molecular distributions. The virtue of the 
simple form Eq. (2) is that it allows the calculation of the reduced density matrix to 
be made explicit. Nevertheless, by perturbative or numerical means, it is presumably 
possible to evaluate the density matrix for the system inside the imaginary sphere for 
an arbitrary density functional S{p}. The result for each eigenmode of the density 
matrix will be formally identical to Eq. (6), although the eigenf unctions \l/ ra will no 
longer be the explicit forms of Eq. (7). Thus, the simple generalized free- field theory 
implied by Eq. (2) is a calculational convenience, rather than an integral part of 
the formalism. Moreover, as is seen below, this simple functional form suffices to 
reproduce not only the necessary area-law form of the geometric entropy (implying 
a constant surface tension), but gives reasonable answers for the surface tension as 
well. 

The next step in the analysis is the calculation of the probability to find the interior 
region in states whose expected density is less than the bulk value. Another point 
must be made clear here. Although the surrounding imaginary sphere represents 
a sharp boundary, the interface itself is not constrained to be as sharp. It is only 
necessary that the boundary of the imaginary sphere be placed outside the interface 
region. Thus, one could imagine making the imaginary sphere extremely large, and 
asking the probability to observe states of density only slightly less than the bulk 
value. For a sufficiently large sphere and small deviation from bulk density, the 
quadratic approximation to the density functional S{p} must therefore be valid. For 
a small sphere, there will however be an error in estimating the size of the bubble of 
low density contained within the sphere. The size of the interface region should be 
dictated by molecular dimensions. (In other words, the existence of an interface in a 
non-critical system should not perturb the system strongly more than a few molecular 
diameters away from that interface). Thus, the ratio of the interface width to the 
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sphere size should vanish for a large enough sphere (which is, of course, the result of 
interest for macroscopic surface tension). Therefore, for a large enough sphere, the 
area of the low-density bubble is asymptotically equal to the area of the surrounding 
imaginary sphere. 

If a different problem is formulated, such as the calculation of the "boundary" 
free energy associated with a bulk liquid in contact with a hard wall, constraints 
on the selected configurations must be introduced. For the boundary free energy 
problem, the density must be zero at the wall. Thus, rather than asking [as in Eq. 
(12)] the probability p an y{N) to find a mode in any state greater than N, one asks 
the probability to find states with a node in the proper location. A more detailed 
analysis would be needed in order to select the proper states to include in this sum. 
However, in any case, there must be factors of p n included in the calculation, and 
thus there should also always be a connection to the geometric entropy. This relation 
should persist even if a better model of a fluid than Eq. (2) is used. 

At least for some problems, it is possible to improve upon the results given by the 
approximate quadratic density functional S{p} by means of a simple perturbative 
analysis, although in general such improvements likely require numerical simulation. 
One example of how a perturbative study might proceed is as follows. Consider a 
more elaborate model of a fluid constructed by replacing the probability functional 
exp[— S^p}] defined by Eq. (2) with a new functional 

exp[-S new {p}] = exp[-S{p}] + e x exp[-S gas {p}], 

where PoS gas {p} — \J[p( r )] 2 dr. Here, e is a small number that gives the relative 
probability of the bulk fluid to be in the gas phase. The new action, S new {p}, has 
two minima, essentially at zero density, p — 0, and at the bulk density, p = po, and 
thus includes a crude characterization of the gas phase of the fluid. One can then 
calculate, using ordinary quantum-mechanical perturbation theory, the change in the 
eigenvalues g n of the density matrix p in due to the additional term in exp[— S{p}} . The 
gas-phase contributions to the geometric entropy —T,g n ln g n can then be expanded 
in a series in the small number e, allowing geometric entropy to be determined in a 
controlled fashion. 

It is therefore clear that no barrier of principle exists in applying the techniques 
developed here to fluid models of arbitrary complexity. The present formalism is thus 
not subject to the easy criticism^ applied to early integral-equation treatments of 
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surface tension. In these early calculations, the direct correlation function appro- 
priate to a bulk fluid was employed essentially to estimate the relative probability 
of various interfacial configurations, whose local density was necessarily far from its 
bulk value. These calculations were thus highly sensitive to the tails of the probabil- 
ity distribution employed, while the bulk direct correlation only provides information 
about the distribution near its peak. The present argument obviates this difficulty 
by counting states, since the specific form of the density matrix for densities far from 
the bulk value contributes little to S in /k = —Tr(pi n \np in ). Thus, by calculating the 
geometric entropy, one is determining the relative probability of the system to be in 
a state whose typical density is outside the region where the probability distribution 
is maximized. Since both the distribution near its peak and the overall normaliza- 
tion of the probability distribution are known, this calculation can be done relatively 
accurately. Moreover, as the density matrix p in is constructed by integrating over all 
fields {p(r)} external to the sphere, correlations of arbitrarily long wavelength are 
included, as they should be! 5 ' 9 ]. Nevertheless, as discussed above, if the effects of 
complex boundary conditions are to be accounted for in a satisfactory way, further 
constraints on the available states must be included. 

4. Evaluation of surface tension 

The calculation of S ou t is easily performed using Eq. (34) of ref. [2]. The 
asymmetric matrix A = M~ 1 N has the same eigenvalues as the symmetric matrix 
A s = M~zNM~z , but is easier to construct numerically. Therefore, the eigenvalues 
of A are computed. By making an expansion in spherical harmonics, the problem 
reduces to finding the eigenvalues A<? of the matrix 

k t (r,r") = - f R T e [l/W}(r,r') T e [WW >") r' 2 dr' (13) 
Jo 

with respect to the measure (r") 2 dr" over the interval {r, r" > R}, where 

i r°° 

T e [W](r',r") = (r'r")-^ kdkW(k) J i+ i(kr') J e+ i(kr") (14) 

Jo 2 2 

and similarly for T^l/W]. Here, W( k) is the Fourier transform of W(r), while 
Je + i (z) is a Bessel function. Then' 2 ! 

OO 1 j j 

Sout/k = E(2^ + l){m(-Aj) + (1 + ln[(l + A,- 1 )^ + A,"]} (15) 
f=o A 
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The sum in Eq. (15) is taken over all eigenvalues Xe for each I. The function 
Poh(k) = 1/W(k) — 1 is taken directly from experiment, obviating the need for an 
independent determination of the density p . As k increases, h(k) goes to zero. Hence, 
the substitutions 

1/W(k) -> p h(k) = w-\k) (16) 

W(k) -> -p ~h(k)/[l + pji{k)\ = w(k) (17) 

were made in Eqs. (13-14) to ensure the convergence of the integrals Eq. (14). Since 
Tg[l](r, r') = (r)~ 2 5(r — r'), subtracting one from W(k) [or 1/TV(A;)] leaves Ag{r, r") 
unaffected for {r, r" > R}, and so the sum Eq. (15) is unchanged by this replacement. 
In contrast to the case^' 3 ! of a massless scalar field theory [where W(k) = k], no 
ultraviolet cutoff is needed here, since h(k) vanishes for large k. Data are generally 
presented as a set of N values of poh(k) at uniformly spaced intervals of k = n x (Ak), 
so the integrals Eq. (14) can be viewed as finite trapezoid- rule sums. The eigenvectors 
of Ag(r,r") are linear combinations of 

u e (k,r") = e(r" - R)(r")-*J t+ i(kr"), (18) 

for values of k contained within the data set, reducing the problem to a numerical 
diagonalization of the N x N matrix 

N 

A e (m,n) = -m[w _1 (mAA;)] ^I e (m,i)[w(iAk)][6 in - il e (i,n)] (19) 

i=i 

for each t. The required indefinite integrals of Bessel functions 

rRAk 

Ie(m,n) = / J e+ i (ms) J e+ i (ns) s ds (20) 
J o 2 2 

can be performed analytically. The computational procedure was checked against the 
test case h(r) = — exp(— mr), for which the integrals Eq. (14) are known exactly. 

The geometric entropy S out (R) for water at 25° C was calculated by applying the 
above procedure to structure function data^ 10 '. The summand of Eq. (15) vanishes 
exponentially with £, albeit with a large decay length (~ 10). By including all values 
of £ < 75, S ou t(R) can be accurately calculated up to R ~ 7. 5 A. Figure 1 displays 
the result. Except for very small R values, S 0Ut (R) is linearly proportional to R 2 , as 
it is for a cutoff scalar field theory^ 2 ' 3 !. Asymptotically, S out (R) ~ 2.50-R 2 , yielding a 
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surface tension of 82 dynes/cm, which is 14% larger than the experimental number 
72 dynes/cm t 11 '. For comparison, scaled particle theory using the usual hard-sphere 
radius of 2.7A predicts' 13 ^ 52 dynes/cm at 303°i^. A detailed molecular dynamics 
simulation t 14 ' of water at 305°fT yields 67 dynes/cm, to be compared with the 71 
dynes/cm obtained experimentally at this temperature. These simulation results 
can be fitted by scaled particle theory, if the significantly larger value of 2.875A 
is chosen for the hard-sphere radius of water. [As there are several implicit scales 
present, dimensional analysis cannot be used to obviate the density-matrix entropy 
calculation. Scaled particle theory, where surface tension is given in terms of the 
hard-sphere radius a by the form 7 = kT / a 2 f (poa 3 ) , is probably the simplest viable 
alternative. Indeed, the only external input to the density matrix calculation, poh(k), 
is dimensionless.] 

The surface tension of a number of other liquids was also calculated using their 
measured structure functions' 15,16 ]. The density matrix results (DM) are given in 
Table 1, along with experimental data I 11 ' 13 ' 17 ! and available estimates' 13 ' from scaled 
particle theory (SPT). The agreement between the density matrix results and exper- 
imental data over a wide temperature range is respectable, especially since the only 
input to the calculation is the (measured) structure factor. By contrast, scaled particle 
theory is quite sensitive to the value chosen for the hard sphere radius. This param- 
eter is not measurable directly, and is usually determined by fitting compressibility 
data. The closest agreement between experiment and the density matrix calculation 
occurs for sodium, potassium, and water, whose structure functions are known with 
the best precision. 

5. Conclusions 

It is worth emphasizing the remarkable fact that the geometric entropy was found 
to be essentially a linear function of the area in every case considered here, in addition 
to the free scalar field theory considered previously by others ' 2 ' 3 '. On very general 
grounds' 3 ' , the eigenvalues of the reduced density matrix (and thus the geometric en- 
tropy) depend upon the area of the boundary surface, and not upon the volume of the 
excluded region. However, it is not clear why the geometric entropy should be linearly 
proportional to this area, especially for situations whose physical basis and mathe- 
matical structure are substantially different. (Thus, for example, the calculation of 
geometric entropy in the "black hole" case required the introduction of an ultraviolet 
cutoff; here, since the structure factor of a fluid vanishes at large wavenumber, no 
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cutoff is needed.) Moreover, the linear dependence of geometric entropy on area pro- 
vides an additional point of support for the present mode of calculation. It not only 
gives the correct "area law" form for this contribution to the free energy, but gives 
reasonable values for the constant of proportionality (which is the surface tension). 

As a rule, few ideas relevant to quantum gravity have practical applications. One 
exception was given here. By determining the geometric entropy of a "black hole" in a 
fluid "field theory" , the liquid surface tension can be evaluated in a novel way. These 
calculations are essentially exact results for a field theory that is an approximation to a 
real liquid. The procedure is therefore automatically self-consistent, which may partly 
account for its plausible accuracy. Possible extensions of the formalism include an 
extraction of the "hydrophobic" potential of mean force, which is a major ingredient 
of protein folding models^ 18 ' . This may be determined by excluding two regions of fluid 
and calculating the resultant geometric entropy as a function of the distance between 
them. The density matrix method is sensitive to the form of the structure factor over 
its entire range. Thus, the method may provide a useful test of approximate theories. 
The major virtue of the density-matrix method is its relative simplicity - molecular 
simulations of entropic quantities are difficult and, hence, rare. Further calculations 
are in progress. 

It is a pleasure to thank E.G.D. Cohen, M.J. Feigenbaum, N.N. Khuri and C.K. 
Zachos for several useful discussions, and Michael Wortis for some particularly tren- 
chant remarks. 
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Table 1. Surface tension 7 (dynes/cm) 





DM 


Experiment 


SPT 


Sodium 


164 (373°K) 


206 (371°K) 


91 (371°K) 


Potassium 


93 (338°K) 


88 (337°K) 


56 (337°K) 


Water 


82 (298°K) 


72 (298°K) 


52 (303°K) 


Chlorine 


31 (298°K) 


18 (293°K) 


20 (293°K) 


Methane 


21 (96°K) 


16 (95°K) 




Nitrogen 


6.4 (77°K) 


9 (77°K) 


9 (77°K) 
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Figure Caption 



Figure 1 

Plot of Sout(R) versus R 2 (A) 2 for water. Dashed line is the function 2.50/? 2 , displayed 
for comparison with plotted data points. 
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